Determining all gas properties in galaxy clusters from the dark 

matter distribution alone 
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ABSTRACT 

We demonstrate that all properties of the hot X-ray emitting gas in galaxy clusters are 
completely determined by the underlying dark matter (DM) structure. Apart from the standard 
conditions of spherical symmetry and hydrostatic equilibrium for the gas, our proof is based on the 
Jeans equation for the DM and two simple relations which have recently emerged from numerical 
simulations: the equality of the gas and DM temperatures, and the almost linear relation between 
the DM velocity anisotropy profile and its density slope. For DM distributions described by the 
NFW or the Sersic profiles, the resulting gas density profile, the gas-to-total-mass ratio profile, 
and the entropy profile are all in good agreement with X-ray observations. All these profiles are 
derived using zero free parameters. Our result allows us to predict the X-ray luminosity profile of 
a cluster in terms of its DM content alone. As a consequence, a new strategy becomes available 
to constrain the DM morphology in galaxy clusters from X-ray observations. Our results can 
also be used as a practical tool for creating initial conditions for realistic cosmological structures 
to be used in numerical simulations. 

Subject headings: dark matter, galaxies: clusters: general, X-rays: galaxies: clusters 



1. Introduction 

Galaxy clusters are the largest equihbrated 
structures in the Universe, consisting mainly of 
dark matter (DM) and hot ionized gas in hydro- 
static equilibrium in the overall potential well. 
Observations of this X-ray emitting gas allow for 
an accurate determination of the properties of 
the dominating DM structure, which can then be 
compared with the results of numerical N-body 
simulations. 

More specifically, the strategy can be outlined 
as follows. In the first place, the equation of hydro- 
static equilibrium can be used to infer the DM den- 



tion (|Navarro et al.||1996| 


Moore et al.||1998 


Die- 


mand et al. 


2004 


Stadel et al. 


2008 


Navarro et 


al. 20081. Moreover, using a very simple connec- 



tion between the gas and DM temperatures which 
has been confirmed by numerical simulations, the 
equation of hydrostatic equilibrium can be com- 
bined with the Jeans equation for the DM to de- 
rive both the DM radial velocity dispersion and 



the velocity anisotropy profile (Hansen & Piffaretti 
2007[|Host et al.||2009[ ). Again, the resulting pro- 
files turn out to be in excellent agreement with 



numerical simulations (Cole & Lacey 1996 Carl- 
berg et al.||1997| ). 



sity profile from X-ray data (Fabricant et al. 1980 1 



Application of this technique in conjunction with 



Measurements of the gas temperature profile 
have demonstrated its virtually universal prop- 



present-day observations 


Voigt & Fabian 


2006 


Pointecouteau et al. 


2005 


yields density profiles 



which are in excellent agreement with those emerg- 
ing from numerical simulations of structure forma- 



erties (1 


)e Grandi & Molendi 


2002 


Kaastra et 


al. 


2004 


Vikhlinin et al. 


2005 


Pointecouteau et 


al. 


2005 


Voigt & Fabian 


2006 


. This universal- 



ity appears surprising since the temperature pro- 
file should encode information about the violent 



1 



gravitational processes taking place during the 
cluster formation, as well as any additional en- 
ergy input, e.g. from a central heat engine, and 
these processes are expected to differ significantly 
from structure to structure. The gas density pro- 
file also exhibits a roughly universal behaviour, 
which has allowed observers to fit remarkably sim- 



ple forms, e.g. a beta-profile (Cavaliere & Fusco- 
Femiano||1976| [S arazm ]|1986| ), to the data. 



A natural question arises as to whether the 
properties of the hot X-ray emitting gas in galaxy 
clusters can be predicted from first principles start- 
ing from the cluster DM distribution alone. 

Attempts in that direction have been pioneered 
by (iMakino et al.||1997| iSuto et al.||1998|). Our 



goal is to go a step further along this direction, 
relaying upon information on the internal cluster 
dynamics which has become available only quite 
recently. 

More specifically, we will show how the gas 
density profile can be obtained directly from the 
underlying DM profile, by combining the equa- 
tion of hydrostatic equilibrium for the gas and 
the Jeans equation for the DM under the stan- 
dard assumption of spherical symmetry. Besides 
the above-mentioned relation between the gas and 
DM temperatures, our derivation rests upon a 
very simple connection between the DM veloc- 
ity anisotropy and the slope of its density profile, 
which has recently emerged in numerical simula- 



tions (Hansen & Moore 2006 Hansen & Stadel 



2006 ) . We thereby demonstrate that the gas den- 



sity profile is completely determined once the grav- 
itationally dominant DM density profile is given. 
Since the gas temperature profile is also known, 
it turns out that the DM distribution dictates all 
the gas properties uniquely. Besides conceptually 
relevant in itself, this fact allows to predict the X- 
ray luminosity profile of a cluster in terms of its 
DM content alone by a method similar to the one 



put forward by Makino et al. (1997); Suto et al. 



(1998 1 . So, a new strategy becomes available to 



constrain the DM morphology in galaxy clusters 
from X-ray observations. Moreover, our findings 
can be employed as a practical tool for creating 
initial conditions for realistic cosmological struc- 
tures to be used in numerical simulations. 



2. Background 

We start by recalling some basic information 
which will be instrumental for our analysis. We 
restrict our attention throughout to regular clus- 
ters, which are supposed to be spherically sym- 
metric and relaxed. The condition of hydrostatic 
equilibrium for the X-ray emitting gas can be writ- 
ten as 

keTg f dlnpg , dlnTg ^ ^ GMtotjr) 

r ' 



/iTTip \ d In r d In 7 



where Pg{r) and Tg{r) arc the gas density and 
temperature profiles, respectively, /i ~ 0.61 is the 
mean molecular weight for the intracluster g as , Trip 
is the proton mass and Aftot('") represents the to- 
tal mass inside radius r. Two conditions have to 
be satisfied in order for Eq. ([ij to hold. First, it 
should be applied to a region considerably larger 
than the gas mean free path, so that local thermo- 
dynamic equilibrium is established. Second, the 
cooling time in that region should be larger than 
the age of the cluster, so that no bulk motion oc- 
curs. The latter condition is generally met outside 
the central region, where the presence of a cooling 
fiow often requires Eq. ([T|) to be replaced by the 
Euler equation (with the velocity term playing a 
nonnegligible role). Because of coUisional relax- 
ation, the gas velocity distribution is isotropic and 
its temperature can be expressed in terms of the 
one-dimensional velocity dispersion cr^ as 



T — 



pmpcrl 



(2) 



Assuming complete spherical symmetry for the 
DM distribution, the two tangential components 
of the DM velocity dispersion, denoted by CTj , are 
necessarily equal, but they are generally allowed to 
differ from the radial component cr^, since DM is 
supposed to be coUisionless. It is usual to quantify 
the DM velocity anisotropy by 



/3 = 1 



^1 



(3) 



and we find it convenient to introduce the mean 
DM one-dimensional velocity dispersion cr^j^ as 



'DM 



3K 



1 



(4) 
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Moreover, in analogy with the case of a gas, we 
also define the DM temperature as [Hansen fc Pif-| 
farettil (pOOTl 



DM 



(5) 



Of course, the coUisionless nature of DM prevents 
any definition of temperature in the thermody- 
namic sense and in fact Tdm is simply meant to 
quantify the average velocity dispersion over the 
three spatial directions. Any completely spher- 
ically symmetric and relaxed DM configuration 
obeys the Jeans equation 

din a1 

dlnr 



dlnpDM 
dlnr 



2/3 



GMtot(r) 



, 
(6) 

where pdmI'') denotes the DM density profile 
( [Binney fc Tremaine|[l987l ). 



3. The temperature profile 

Early studies of the X-ray emission from regu- 
lar clusters were based on the assumption of an 
isothermal gas distribution, simply because the 
Einstein observatory and ROSAT were unable to 
determine the cluster temperature profiles. The 
observed X-ray emission is produced by thermal 
bremsstrahlung ( Sarazin|[l986 l, so for Tg — const, 
it follows that Pg(r) is proportional to the square 
root of the deprojected X-ray surface brightness. 
In such a situation, a good fit to the data was 



provided by the beta-model (Cavaliere & Fusco- 
fFcmianO||1976| |Sarazin||1986[ ) 



Pg(0) 







1 + 









3Ait/2 



(7) 



where ax < 0.5 Mpc is the X-ray core radius. Note 
that /?fit has nothing to do with the DM velocity 
anisotropy. Typically, most of the emission comes 
from the region r > 0.5 Mpc, and so Eq. (|7| can 
be approximated by the power-law 

pg{r)c,pg{0){ — ] 



Now, by inserting Eq. (|8| and Tg 
Eq. Q, we find 



(8) 



const, into 



Mtot(r) = 



G 



(9) 



where Eq. ([2| has been used. As is well known, un- 
der the assumption of isotropic velocity distribu- 
tion (/3 = 0), a mass profile of the form M(r) esc r 
describes a singular isothermal sphere (SIS) model 
in which the velocity dispersion is everywhere con- 



stant (Binney & Tremaine 19871. Denoting by a 



the one-dimensional velocity dispersion, we explic- 
itly have M(r) = {2(j'^/G) r. Owing to the fact 
that the leading contribution to Aftot(^) comes 
from DM, it follows that Mtot(r) ~ Muuir)- As 
a consequence, Eq. ([9| can be rewritten as 



(10) 



and the comparison of Eqs. ([9| and ( 10 1 entails in 
turn 



'DM 



1.5/3fit erf, 



(11) 



Observations performed with the Einstein obser- 
vatory and ROSA T yield 0.5 < /3fit < 0.9 with a 
median /3fit ^ 0.67 ( [Bahcall fc Lubin|1994[ ). Thus, 
on average we get 



2 

"DM 



which implies 



Tt 



DM 



g ' 



(12) 



(13) 



thanks to Eqs. ^ and 

Only with the advent of the ASCA and Beppo- 
SAX satellites did it become possible to measure 
the cluster temperature profiles, which turned out 
to be described by a polytropic gas distribution 
to first approximation. Higher-quality data are 
currently provided by Chandra and XMM-Newton 
satellites, which have shown that the gas temper- 
ature profiles possess a very simple and nearly 



universal behaviour (see Vikhlinin et al.| (20061 



for a thorough discussion). Basically, it increases 
rapidly from a small (possibly non-zero) value 
in the centre, to a maximum at a radius about 
O.lriso, and then declines slowly by a factor of 
2 — 3 at (0.6 — 0.8) rigQ. Here, rigo is defined as 
the radius within which the mean total density is 
180 times the critical density at the redshift of the 
cluster. The necessary X-ray background subtrac- 
tion makes it very difficult to accurately measure 
the temperature further out. 

As mentioned above, our main goal is the de- 
termination of the gas density profile Pg{r) once a 
specific dark matter distribution AfDM('') is given. 



3 



Supposing as before that Mtot{r) — Mj^uir), it is 
evident that Pg{r) follows from Eq. ([T| provided 



that TJr)js specified. Previous studies (IMakino 



et al.| 1997 |Suto et aL 1998 1 accomplished this 
task by assuming 



(14) 



which was suggested to formalize the condition 
that the gas temperature is close to the virial tem- 
perature of the DM. However, the virial theorem 
is a global relation that characterizes a cluster as 
a whole - it just arises by integrating the Jeans 
equation over the system - and so its validity is 
open to criticism. Of course, in the lack of any 
other information about the gas temperature this 
was the only viable possibility. However, recent 
progress allows for a considerable improvement on 
this point. 

As a matter of fact , this stumbling block can be 
side-stepped in a remarkably simple fashion. Be- 
cause of the equivalence principle, the velocity of a 
test particle in an external gravitational field is in- 
dependent of the particle mass. This circumstance 
leads to the guess 



(15) 



This relation was tested against numerical sim- 
ulations (Host et al. 2009), which demonstrated 



its validity with k = 1 to a very good approxi- 



mation. These numerical simulations (Kay et al. 
2007[ [Springcl .200"5| |Valdarnini. 200'6| are reliable 
only on scales greater than ~ 0.1r25oo, while the 
best X-ray observations are sensitive to a radius 
which is almost a factor 3 smaller. It is there- 
fore possible that heating or cooling may shift k 
away from unity in the very centre. Hence, out- 
side that region k = 1 is expected. Actually, a look 



back at Eq. ( 13 ) confirms the remarkable fact that 
K = 1 holds regardless of the actual shape of the 
DM velocity anisotropy profile (3{r). As we shall 
see, starting from a specific underlying DM density 
profile puuir), one can evaluate T^uir) and then 
get the gas temperature profile Tg{r) uniquely. 

Before closing this section, a remark is in order. 
Observations show that some clusters lack a cen- 
tral cooling flow. In such a situation, hydrostatic 
equilibrium is expected to hold all the way down 
to the centre. Actually, for typical central values 



of the electron number density rig ~ 1cm ^ and 
temperature T ~ 10*K ~ 8.5keV ( |Sarazin|[l986| , 
the scattering time turns out to be tgcat ^ lO^yr, 
which is much smaller than the corresponding gas 
cooling time tcooi ~ 10''yr, so that local hydro- 
static equilibrium is indeed fulfilled outside a cen- 
tral spherical region of radius ~ 1 pc. Assuming 
further that the gas temperature is roughly con- 
stant in the inner cluster region, the gas density 
profile cannot be cuspy as long as Mdm(?') oc r" 
with a > 1 for r ^ 0. This is at odds with blind 
extrapolations of fitting formulae for the temper- 
ature and density such as those used in |Vikhlinin| 
eTaLl ((20061 . 



4. The density profile 

We now proceed to the actual derivation of the 
gas density profile Pg{r) from the properties of the 
dominating DM distribution. 

As a preliminary step, we notice that Eqs. ([T]) 
and (pi can be trivially combined to yield 



ksTg / d In pg ^ d In Tg ^ 
fjinip \ d In r d In r y 
^ d In pdm d In crj? 



d In r d In r 



-2/3 . (16) 



Owing to Eqs. (|5| and (15 1 with k = 1, straight- 
forward manipulations permit to recast Eq. (16 1 
into the form 



7DM + 2/?-|- |/3- 



2 d/3 



din erf 



d In 7' 3 d In r 



(17) 



where we have defined the density slopes 7dm('') 
of the DM and 7g(r) of the gas as 



lx{r) 



din px 
dlnr 



(18) 



with X standing for either DM or g. We stress that 



Eq. (17 1 captures a crucial point of the present 



investigation: only the gas density slope appears 
on its left-hand side, whereas only quantities per- 
taining to the DM appear on its right-hand side. 
It should be appreciated that this result merely 
relies upon the equality of gas and DM tempera- 
tures and - unlike in previous studies ( Cavaliere k. 



Fusco-Femiano||1976| [Makino et ar]|1997| |Suto et 

al.|1998 ) - no assumption is being made about the 



actual gas temperature structure (e.g. isothermal 
or polytropic). 
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Fig. 1. — The derived gas density profile, as- 
suming that Pg/pBM = 10% at ro, which is the 
scale length of the NFW profile. The upper curve 
(black) is the DM density, and the 3 lower lines 
show gas profiles modelled with extreme varia- 
tions in the possible DM velocity anisotropy (green 
dot-dashed is isotropic (/3 = 0), red sohd is using 
/? = - 0.2 (7 0.8) ( [Hansen fc Stadell[2006l) , an d 
blue dashed is using j3 = -0.13 7 ( |Hansen|2008| ). 



Next, we use the fact that the DM anisotropy 
profile f3{r) turns out to be almost linearly related 
to the slope of the DM density profile 7dm (^)- 
This result has been obtained from numerical sim- 
ulations and holds with a scatter of about 0.05 



(Hansen & Moore 2006 Hansen & Stadel 20061. 



It has recently been confirmed by high-resolution 
numerical simulations (Navarro et al. 2008) and 



moreover it has been derived analytically ( Hansen 



2008) (see also Zait et al. (2008); Wojtak et al 



(2008);|Sah^dor-Sole et aL (2007)) 



Getting the gas density profile Pg{r) now in- 
volves a few simple steps. Our only input is the 
DM density profile puuif), like e.g. an NFW pro- 
file. Thanks to Eq. ( [l8| , we rewrite the Jeans 
equation (|6| as 



+ <j'r (7DM(r) + 2/3(r)) + Illll^ ^ Q , 

(19) 

whose solution is easily found to be 



GAftot(r) 



G 

W) 



'^,Biy)M^^ (20) 



Fig. 2. — The derived slope of the gas density pro- 
file, assuming an NFW profile for the DM. Same 
notation as in Figure [T| Th e inner point is taken 



from Vikhlinin et al. (2006) and the three outer 



points are taken from Ettori & Balestra ( 2008 1 



ith 

i?(r) = pDM(r)exp(- 2 /"dr'^ (21) 



Using the relation between f3{r) and 7dm('"), we 



finally obtain the gas density profile from Eqs. ( 17 1 



and (18) 



In practice, such a procedure can be imple- 
mented iteratively. In first approximation, we as- 
sume that the gas mass contribution is negligible, 
so that we have Mtot(f) = Mdm(?')- lii the next 
iterations, we include the gas mass in the calcula- 
tion of (7^ (r) . Although the gas mass is taken into 
account perturbatively, any desired accuracy can 
be achieved by a sufficient number of iterations. 

An example of the application of this strategy 
is shown in Figure [l] where the DM density is as- 
sumed to follow an NFW profile (black solid line) . 
The three lower lines are the gas density profiles 
obtained with a range of different possible DM ve- 
locity anisotropy profiles. The details of the gas 
density profile are easier seen in the slope, which 
is shown in Figure |2] Note that for an inner DM 
slope of about 1 (in agreement with the observa- 
tions (Voigt & Fabian 2006 Pointecouteau et al 



2005)) the inner gas slope should also be close to 
1. This is in good agreement with the fits from 
Vikhlinin et al. ( j2006) , which have an average of 
0.8 for the extrapolated inner slope. Also the 



5 



slopes found by Ettori & Balestra (20081 agrees 
with an NFW profile. 




0L_ 
0.01 



0.10 1.00 
Radius (r / Tq) 



Fig. 3. — The slope of the gas density profile, as- 
suming a Sersic profile with n — 5 for the DM. 
Same notation as in Figures [T] and 2. 

A widely used alternative to the NFW profile 
is the Sersic (or Einasto) profile, which generalizes 
the de Vaucouleurs profile traditionally used to fit 
the optical surface brightness of elliptical galaxies. 
It has been shown that the Sersic profile models 
the deprojected DM density at least as well as the 



NFW ( 


Navarro et al. 


2004 


Merritt et al. 


2006 


Salvador-Sole et al. 


2007 


). This profile contains 



3 free parameters: two scaling constants for the 
density and the radius - po ^tnd respectively - 
and one shape parameter n 



P{r) = Po exp 




(22) 



The constant 6„ is a function of the index n and is 



tabulated e.g. by Mazure fc Capelato (2002). The 



radial velocity dispersion derived from the Ser- 
sic profile has, like the NFW profile, the property 
of reaching its maximum near rg where the slope 
is 7 = 2. In Figure |3] we present the gas density 
slope, assuming a Sersic profile for the underly- 
ing DM density. There is not sufficient statistical 
power in the data to discriminate between the un- 
derlying DM density and/or velocity anisotropy 
profiles from this analysis. 

Both the NFW and the Sersic profile are con- 
sistent with observations because they have the 
appropriate slope in the inner and outer observed 



region. Since we cannot exclude one or the other 
by relying upon their shape, we choose the NFW 
model for the underlying DM in the rest of our 
treatment. 

Since the gas density profile differs from the un- 
derlying DM density profile, there will also be a 
radial variation in the local and cumulative gas 
fractions, which are defined as 



and 



Ptot(r-) 



(23) 



(24) 



respectively. In order to test this in more detail, 
we used the 16 clusters analysed in [Host et aL] 
( 2009 ) , which is a selection of highly relaxed clus- 



ters at both low and intermediate redshifts (Kaas 



tra et al.|[2004l [Piffaretti eral][2005l [Morandi et 
al.||2007[ ) observed with XMM-Newton and Chan- 
dra. Under the assumption of hydrostatic equi- 
librium, we find the local gas fraction exhibited 
in Figure |4j The local gas fraction clearly in- 
creases as a function of radius, which demonstrates 
that the DM velocity anisotropy cannot vanish 
(green dot-dashed line in Figure |4|. The gas den- 
sity fraction roughly increases as a power-law in 
radius and we have approximately fg{r) r^'^. 
The solid (red) and dashed (blue) lines are for 
NFW DM profiles, with different radial DM ve- 
locity anisotropics. From Figure |4] there is a clear 
difference between the data and the predictions 
in the outer region, which may be due to an un- 
derestimation of the total mass due to breakdown 



of hydrostatic equilibrium (Piffaretti & Valdarnini 



2008) 



It is important to keep in mind that these pro- 
files do not contain any free parameters. Every 
quantity is calculated from the dark matter distri- 
bution alone. 

The agreement in the inner region is better and 
it should be kept in mind that different DM den- 
sity and velocity anisotropy profiles give rise to 
different curves. It may therefore be possible to 
use the shape of 4>g{r) to recover these DM pro- 
files in the future. 

As a further step, we discuss some of the im- 
plications of our main result. Indeed, with a full 
description of the gas that is directly derived from 
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Fig. 4. — The observed (f>g from 16 relaxed galaxy 
clusters. We have scaled the gas mass fraction 
with free parameters, to make the radial variation 
more visible. There is a rough trend that the gas 
mass fraction increases as r°'^, which is in clear 
disagreement with the assumption that the DM 
velocity anisotropy should vanish. The 3 curves 
all assume an NFW density profile, and different 
assumed connections between the DM slope and 
velocity anisotropy, as in Figure [l] 

the dark matter distribution, we can predict addi- 
tional observable quantities besides the gas frac- 
tion described above. 

One of these quantities is the entropy, which is 
often characterized by the adiabatic coefficient Kg 
of the gas 



/ 



-2/3 



(25) 



Our previous results entail that these profiles are 
almost perfect power laws regardless of the P pro- 
file. The slope changes slightly for the different 
/? profiles (between 1.1 and 1.3). This theoretical 
prediction is in good agreement with the X-ray 
observations, which generally produce power-law 



entropy profiles ( 


Piffaretti et al. 


2005 


Pratt, G. 


W. & Arnaud, M.||2005 Donahue et al. 


20061. 



Another quantity that we are able to predict is 
the gas X-ray emissivity e. As a matter of fact, 
e can be estimated either analytically - because 
e oc n^T 2 - or by numerical codes like MeKaL 



(Mewe et al. 19851 in order to include the line 
emission contribution. The latter strategy is es- 
pecially well suited for cooler clusters, because a 



substantial amount of their luminosity stems from 
emission lines. On the other hand, the luminosity 
of hotter clusters is dominated by the continuum 
emission. In either case, the surface brightness 
can be inferred from a given DM profile and this 
can in turn be compared with observations. In 
this way, it is possible to construct an algorithm 
that adjusts the proposed DM profile until the sur- 
face brightness best-fits observations and thereby 
single out the optimal DM profile. It is clear 
that this method will be limited due to dumpi- 
ness, which may be crucial in the outer parts, and 
one may have to combine the X-ray analysis with 
SZ observations to correctly extract all the pro- 
files. We hope to discuss these aspects in more 
detail in the future. 

Whereas numerical simulations have demon- 
strated that the gas and dark matter temperatures 
are equal in large parts of a galaxy cluster, they 
cannot probe the very centre of the clusters. It 



is therefore possible that n in Eq. (15) departs 
from unity as r if there is significant cooling 
or heating. However, in our derivation of the gas 
density profile we have assumed k = I everywhere. 

It goes without saying that we can turn the 
argument around and use the observed gas profile 
to determine k. Basically, we can insert Eq. ( 15 1 



into Eq. (16 1 and solve for n. Furthermore, since 
we are now interested in the central region where 
(3 is likely to be vanishingly small, we discard all 
terms involving P in the resulting equation. So, in 



place of Eq. (17 1 we presently get 
7g -I- d In /d In r 



7dm -I- d In In r 



(26) 



which in principle allows to measure k directly 
from X-ray observations. Such measurement can 
be used or tested in future numerical simulations 
when the increased particle number will allow sim- 
ulations to probe closer to the cluster centre. 

5. Conclusions 

We have shown that all properties of the hot 
X-ray emitting gas in galaxy clusters are com- 
pletely determined by its underlying DM struc- 
ture. Apart from the standard conditions of 
spherical symmetry and hydrostatic equilibrium 
for the gas, our proof is based on the Jeans equa- 
tion for the DM and two simple relations which 
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have recently emerged from numerical simulations. 
One is the equality of the gas and DM tempera- 
tures. The other is an almost linear relationship 
between the DM velocity anisotropy profile and its 
density slope. For DM distributions described by 
the NFW or the Sersic profiles, the resulting gas 
density profile, the gas-to-total-mass ratio profile 
and the entropy profile are all in good agreement 
with X-ray observations. We feel that our result is 
conceptually relevant in itself. Moreover, it allows 
to predict the X-ray luminosity profile of a clus- 
ter in terms of its DM content alone. Therefore, 
a new strategy becomes available to constrain the 
DM morphology in galaxy clusters from X-ray ob- 
servations ( Frederiksen et al.|2009 ). This strategy 
may constrain morphology parameters better be- 
cause of the tighter errors on surface brightness, 
but requires the structure to be very relaxed, and 
ignores dumpiness, and thus cannot be used on 
every cluster. Our results can also be used as 
a practical tool for creating initial conditions for 
realistic cosmological structures to be used in nu- 
merical simulations. 
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